First Measurement of the p Spectral Function in High-Energy Nuclear Collisions 
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We report on a precision measurement of low-mass muon pairs in 158 AGeV indium-indium 
collisions at the CERN SPS. A significant excess of pairs is observed above the yield expected 
from neutral meson decays. The unprecedented sample size of 360 000 dimuons and the good mass 
resolution of about 2% allow us to isolate the excess by subtraction of the decay sources. The shape 
of the resulting mass spectrum is consistent with a dominant contribution from tt + -k~ — » p — > p + p~ 
annihilation. The associated space-time averaged p spectral function shows a strong broadening, 
but essentially no shift in mass. This may rule out theoretical models linking hadron masses directly 
to the chiral condensate. 



PACS numbers: 25.75.-q, 12.38. Mh, 13.85.Qk 



According to quantum chromodynamics (QCD), 
strongly interacting matter will, at sufficiently high tem- 
peratures or baryon densities, undergo a phase transi- 
tion from a state of hadronic constituents to quark mat- 
ter, a plasma of deconfincd quarks and gluons. At the 
same time, chiral symmetry, spontaneously broken in the 
hadronic world, will be restored. High-energy nucleus- 
nucleus collisions provide the only way to investigate this 
issue in the laboratory. 

Among the observables used for the diagnostics of the 
hot and dense fireball formed in these collisions, lepton 
pairs are particularly attractive. In contrast to hadrons, 
they directly probe the entire space-time evolution of 
the fireball and freely escape from the interaction zone, 
undisturbed by final-state interactions. In the hadronic 
phase, thermal dilcpton production in the mass region <1 
GeV/c 2 is largely mediated by the light vector mesons p, 
ui and </>. Among these, the p(770 MeV/c 2 ) is the most 
important, due to its strong coupling to the 7T7t chan- 
nel and its lifetime of only 1.3 fm/c, much shorter than 
the lifetime of the fireball. These properties have given 
it a key role as the test particle for "in-medium modi- 
fications" of hadron properties close to the QCD phase 



boundary. The "dual" role of the p, probing the medium 
and having the medium probe its own properties, was 
pointed out already more than 20 years ago: changes 
both in width and in mass were suggested as precursor- 
signatures of the chiral transition [l| . There seems to be 
some consensus now that the width of the p should in- 
crease towards the transition region, based on a number 
of quite different theoretical approaches 
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the other hand, no consensus exists on how the mass of 
the p should change in approachingthe transition: pre- 
dictions exist for a decrease 0IHJJI1], a constant behav- 
ior 0, 0] , and even an increase . Strongly increasing 
widths and nearly constant masses are also predicted for 
pions and for nucleons 0- 

Experimentally, low-mass electron pair production 
was previously investigated at the CERN SPS by the 
CERES/NA45 experiment for p-Be/Au, S-Au and Pb- 
Au collisions [ToL 111! H^ . The common feature of all 
results from nuclear collisions was an excess of the ob- 
served dilcpton yield above the expected electromagnetic 
decays of neutral mesons, by a factor of 2-3, for masses 
above 0.2 GeV/c 2 . The surplus yield has generally been 
attributed to direct thermal radiation from the fireball, 
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dominated by pion annihilation 7r + 7r~ — > p — > 
with an intermediate p which is strongly modified by the 
medium. Statistical accuracy and mass resolution of the 
data were, however, not sufficient to reach any sensitivity 
to the character of the in-medium changes. 

In this Letter, we present first results on low-mass 
muon pair production from the NA60 experiment at the 
CERN SPS. Compared to the results from CERES, we 
have been able to improve the statistical accuracy by a 
factor of >1000 and the mass resolution by a factor of 2- 
3. This has allowed us, for the first time in this field, to 
isolate the excess by subtraction of the expected sources. 
The essential features of the resulting mass spectrum, in- 
terpreted as the space-time averaged p spectral function 
associated with pion annihilation, can be summarized in 
a few words: the p strongly broadens, but does not show 
any shift in mass. 

The NA60 apparatus complements the muon spec- 
trometer previously used by NA50 with a high- 
granularity silicon pixel telescope of unprecedented ra- 
diation tolerance 0, The telescope, embedded 
in a 2.5 T dipole magnet in the vertex region, tracks 
all charged particles before the hadron absorber and 
determines their momenta independently of the muon 
spectrometer. The matching of the muon tracks be- 
fore and after the hadron absorber, both in angular 
and momentum space, improves the dimuon mass res- 
olution in the region of the vector mesons u>, <fi from ^80 
to ~20 MeV/c 2 , significantly reduces the combinatorial 
background due to it and K decays and makes it possible 
to measure the muon offset with respect to the inter- 
action vertex [J^. The additional bend by the dipole 
field in the target region deflects soft muons into the ac- 
ceptance of the muon spectrometer, thereby strongly en- 
hancing the opposite-sign dimuon acceptance at low mass 
and low transverse momentum with respect to all previ- 
ous dimuon experiments [Hi ]. The rapidity coverage is 
3.3<y<4.3 for the p, at low p-r, compared to 3<y<4 for 
the J/ip. The acceptance, about 1% or smaller relative 
to 4-7T in the whole low- mass region 0], increases with 
Pt by a factor of 3 from to 2 GeV/c for the p and, inte- 
grated over pr, increases with mass by factors of 3-6 from 
0.4 to 0.8 GeV/c 2 , depending on the input p-r spectrum. 
Finally, the selective dimuon trigger and the fast readout 
speed of the vertex tracker allow the experiment to run 
at very high luminosities, leading to an unprecedented 
level of statistics for low-mass lepton pairs. 

The results reported in this Letter were obtained from 
the analysis of data taken in 2003 with a 158 AGcV 
indium beam, incident on a segmented indium target 
of seven disks with a total of 18% (In- In) interaction 
length. At an average beam intensity of 5T0 7 ions per 
5 s burst, about 3-10 12 ions were delivered to the ex- 
periment. A total of 230 million dimuon triggers were 
recorded on tape, collected with two field settings of the 
muon-spectromcter magnet; the data reported here cor- 



respond to the lower of the two settings. 

The data reconstruction starts with the muon- 
spectrometer tracks. Next, pattern recognition and 
tracking in the vertex telescope are done; the interaction 
vertex in the target is reconstructed with a resolution 
of ~200 p,m for the z-coordinate and 10-20 p,m in the 
transverse plane. Only events with one vertex arc kept; 
interaction pileup and reinteractions of secondaries and 
fragments are rejected. Finally, each muon-spectromcter 
track is extrapolated to the vertex region and matched 
to the tracks from the vertex telescope. The matching 
is done using the square of a weighted distance (x 2 )> 
in a space of angles and inverse momenta, between the 
two tracks 0. Correct matches are statistically dis- 
tinguished from fake matches (associations of muons to 
non-muon vertex tracks) by the very different shape of 
the matching x 2 distributions. Matched tracks are kept 
if their x 2 is below a certain threshold. Pair reconstruc- 
tion efficiencies vary by <10% both vs. mass and vs. 
centrality. 

The combinatorial background of uncorrelatcd muon 
pairs mainly originating from n and K decays is deter- 
mined using a mixed- event technique jl5l |. Two single 
muons from different like-sign dimuon triggers are com- 
bined into muon pairs in such a way as to accurately ac- 
count for details of the acceptance and trigger conditions. 
The quality of the mixed-event technique can be judged 
by comparing the like-sign distributions generated from 
mixed events with the measured like-sign distributions. 
It is remarkable that the two agree to within ~1% over 
a dynamic range of 4 orders of magnitude in the steeply 
falling mass spectrum jl5l | . After subtraction of the com- 
binatorial background, the remaining opposite-sign pairs 
still contain "signal" fake matches, a contribution which 
is only 7% of the combinatorial background level. It has 
been determined in the present analysis by an overlay 
Monte Carlo method. We have verified that an event- 
mixing technique gives the same results, both in shape 
and in yield, within better than 5%. More details on the 
experimental apparatus and data analysis will be given 
in a forthcoming extended paper; for now see [l5L Uj\ 

Fig. n shows the opposite-sign, background and signal 
dimuon mass spectra, integrated over all collision cen- 
tralizes. After subtracting the combinatorial background 
and the signal fake matches, the resulting net spectrum 
contains about 360 000 muon pairs in the mass range of 
the figure, roughly 50% of the total available statistics. 
The average charged-particle multiplicity density mea- 
sured by the vertex tracker is dN c h/dn =120, the average 
signal-to-background ratio is 1/7. For the first time in 
nuclear collisions, the vector mesons uj and <j) are com- 
pletely resolved in the dilepton channel; even the rj—>pp 
decay is seen. The mass resolution at the oj is 20 MeV/c 2 . 
The subsequent analysis is done in four classes of colli- 
sion centrality defined through the charged-particle mul- 
tiplicity density: peripheral (4-30), scmipcripheral (30- 
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FIG. 1: Mass spectra of the opposite-sign dimuons (up- 
per histogram), combinatorial background (dashed), signal 
fake matches (dashed-dotted), and resulting signal (histogram 
with error bars). 

110), semicentral (110-170) and central (170-240). The 
signal-to-background ratios associated with the individ- 
ual classes are 2, 1/3, 1/8 and 1/11, respectively. 

The peripheral data can be described by the expected 
electromagnetic decays of the neutral mesons 0J. Muon 
pairs from resonance (77, p, u>, <j>) and Dalitz (77,77 de- 
cays were simulated using the improved hadron decay 
generator GENESIS as the input and GEANT for 
transport through the detectors. Four free parame- 
ters (apart from the overall normalization) were used 
in the fit of this "hadron decay cocktail" to the periph- 
eral data: the cross section ratios rj/u), p/uo and <j>/ui, 
and the level of D meson pa ir decays; the ratio 77 / 77 
was kept fixed at 0.12 0, 0|. The fits were indepen- 
dently done in 3 bins of dimuon transverse momentum: 
Pt<0.5, 0.5<pt<1 and pr>l GeV/c. The particle ra- 
tios obtained from the fits without p T selection, corrected 
for acceptance and extrapolated to full phase space are 
t ? /cj=0.88±0.03, 0/w=O.O94±O.OO4 and p/w=1.62±0.10. 
The errors quoted are purely statistical; the systematic 
errors are estimated to be of order 10%, dominated by un- 
certainties of the branching ratios (using the much more 
accurate value for u> — > e + e~ instead of p + p~ with the 
assumption of e/j, universality). The value for 77/a; agrees, 
within err ors, with the literature average of 0.82±0.11 for 
p-p, p-Be [l(j|. The value for <f>/uj is subject to some <p- 
enhancement already in peripheral nuclear collisions and 
can therefore not be directly compared to p-p values. 
Within 10% these two ratios do not depend on the pair 
Pt- In contrast, the particle ratio p/tu monotonically de- 
creases with pt reaching 1.20±0.09 at pt>1 GeV/c, 
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FIG. 2: Isolation of an excess above the electromagnetic de- 
cays of a "conservative" hadron decay cocktail (see text). To- 
tal data (open circles), individual cocktail sources (solid), dif- 
ference data (thick triangles), sum of cocktail sources and 
difference data (dashed). 



where it a gree s with the literature average of 1.0±0.15 for 
p-p, p-Be [ljj within errors. This suggests that some 7T7T 
annihilation, enhancing the yield of the low-pr p, con- 
tributes already in peripheral collisions (see below) . The 
level of charm decays obtained from the fits is essentially 
determined by the measured yield in the mass interval 
1.2<M<1.4 GeV/c 2 . The results altogether, and in par- 
ticular those on the pr independence of n/uj and <j)/u>, 
indicate that the acceptance of NA60 both in mass and 
Pt, including the critical low-mass (77-Dalitz), low-p^ re- 
gion (<0.5 GeV/c), is reasonably well understood. Apart 
from the sources discussed, no further input was required 
to describe the peripheral data. 

In the more central bins, the data can no longer be de- 
scribed on the basis of the standard hadron decay cocktail 
alone, but is indicative of the existence of an excess yield. 
Since the particle ratios are expected to be different from 
the peripheral data, global fits to the more central data 
are bound to bias both the extracted cocktail parameters 
and an excess with a priori unknown characteristics. We 
have therefore used a novel procedure, made possible by 
the high data quality; it is illustrated in Fig. The ex- 
cess is isolated by subtracting the cocktail, without the 
p, from the data. The cocktail is fixed, separately for the 
major sources and in each centrality bin, by a "conserva- 
tive" approach. The yields of the narrow vector mesons 
cu and <j> are fixed so as to get, after subtraction, a smooth 
underlying continuum. For the 77, an upper limit is de- 
fined by "saturating" the measured data in the region 
close to 0.2 GeV/c 2 ; this implies the excess to vanish at 
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FIG. 3: Excess mass spectra of dimuons. The cocktail p 
(solid) and the level of uncorrelated charm decays (dashed) 
are shown for comparison. The errors are purely statistical; 
see text for systematic errors. 



12%, 25% and 25% in the 4 centrality bins, from periph- 
eral to central. Uncertainties associated with the hadron 
decay cocktail and its subtraction including branching 
ratios, transition form factors and helicity distributions 
have been estimated to be around 15%. Since the back- 
ground causing the dominant errors is essentially flat, 
the p-like structure above the continuum is much more 
robust. 

The qualitative features of the mass spectra in Fig. [3] 
are consistent with an interpretation of the excess as 
dominantly due to tttt annihilation. Among the many 
different theoretical predictions for the properties of the 
intermediate p mentioned in the introduction, only two 
have been brought to a level suitable for a quantitative 
comparison to the data: the broadening scenario of 0,0 
and the moving- mass scenario related to [(J 0] • Both are 
evaluated for In-In at dN c h/drj=lAQ within the same fire- 
ball evolution, taking explicit account of temperature as 
well as of baryon density [2£j. In Fig. 01 these predic- 



very low mass, by construction. The r\ resonance and 
lu Dalitz decays are now bound as well; r\ jr\ is fixed as 
before. The cocktail p (only required in Figs. [3] and 0] 
for illustration purposes) is bound by the ratio p/ui=l.2, 
found at high y>t (>l-6 GeV/c) for all centralities. The 
accuracy in the determination of the u and <p yields by 
this subtraction procedure is on the level of 1-2%, due 
to the remarkable local sensitivity, and not much worse 
for the rj. The qualitative features of the resulting differ- 
ence spectrum are robust towards yield changes even on 
the level of 10%, again because the consequences of such 
changes are highly localized. 

The excess mass spectra for all 4 multiplicity bins, re- 
sulting from subtraction of the "conservative" hadron 
decay cocktail from the measured data, are shown in 
Fig. The cocktail p and the level of charm decays, 
found in the 3 upper centrality bins to be about 1/3 
of the measured yield in the mass interval 1.2<M<1.4 
GeV/c 2 E3, are shown for comparison. The qualitative 
features of the spectra are striking: a peaked structure 
is seen in all cases, broadening strongly with centrality, 
but remaining essentially centered around the position 
of the nominal p pole. At the same time, the total yield 
increases relative to the cocktail p. their ratio reaching 
values above 4 for M<0.9 GeV/c 2 in the most central 
bin. Such values are consistent with the results found by 
CERES ^3] , if the latter arc also referred to the cock- 
tail p and rescaled according to the different multiplicity 
density. The errors shown are purely statistical. The 
dominant sources of systematic errors are connected to 
the uncertainties in the levels of the combinatorial back- 
ground (1%) and fake matches (5%). On the basis of 
these values and the signal-to-background ratios, the sys- 
tematic errors in the broad continuum region 0.4<M<0.6 
and 0.8<M<1.0 GeV/c 2 are estimated to be about 3%, 
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FIG. 4: Comparison of the excess mass spectrum for the 
semicentral bin to model predictions, made for In-In at 
dNchl 'dr?=140. Cocktail p (thin solid) , unmodified p (dashed), 
in-medium broadening p 0,0 (thick solid), in-medium mov- 
ing p related to H0 (dashed-dotted). The errors are purely 
statistical. The systematic errors of the continuum are about 
25% (see text). The open data points show the difference 
spectrum resulting from a decrease of the rj yield by 10% 
(which should also be viewed as a systematic error). 

tions (as well as the unmodified p) arc confronted with 
the data for the semicentral bin (same charged-particlc 
rapidity density). Note that the integrals of the theoreti- 
cal spectra are normalized to the data in the mass interval 
M<0.9 GeV/c 2 . The unmodified p is clearly ruled out. 
The specific moving-mass scenario plotted here, which 
fitted the CERES data 0,0], is also ruled out, showing 
the much improved discrimination power of the present 
data. The broadening scenario appears more realistic. 
However, the nearly symmetrical broadening around the 
p pole seen in the data is not reproduced by this model. 
The remaining excess at M>0.9 GeV/c 2 may well be re- 
lated to the prompt dimuon excess found by NA60 in 
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the intermediate mass region [15|. Processes other than 
2ir, i.e. Air ... could possibly account for the region 
M>0.9 GeV/c 2 

The data shown in this Letter have not been corrected 
for the mass- and py-dependent acceptance of the NA60 
setup. The theoretical calculations shown in Fig. 0] were 
therefore also propagated through the acceptance filter to 
allow for fair comparison with the data. It is interesting 
to note that, by coincidence, the propagation of theoret- 
ical calculations based on a white spectral function (e.g. 
qq annihilation |20l| ) yields a mass spectrum flat within 
10% up to about 1 GcV/c 2 , without any bump struc- 
ture. In other words, the always existing steep rise of 
the theoretical input at low masses, due to thephoton 
propagator and a Boltzmann-like factor 0, IE EI 01 ; ^ s 
just about compensated by the falling acceptance in this 
region as long as no pr cut is applied. The data and 
model predictions shown in Figs. and 0] can therefore 
be interpreted as spectral functions of the p, averaged 
over momenta and the complete space-time evolution of 
the fireball. The flat part of the measured spectra may 
thus reflect the early history close to the QCD boundary 
with a nearly divergent width, while the narrow peak on 
top may just be due to the late part close to the thermal 
freeze-out, approaching the nominal width. 

We conclude, in quite general terms and independently 
of any comparison to theoretical modeling, that the p 
primarily broadens in In-In collisions, but does not show 
any noticeable shift in mass. More detailed work includ- 
ing precise p T dependences is under way to consolidate 
these findings. 
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